Andreev drag effect in ferromagnetic-normal-superconducting systems 
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We investigate conductances and current correlations in a system consisting of a normal mul- 
tichannel conductor connected to one superconducting and two ferromagnetic electrodes. For an- 
tiparallel orientation of the ferromagnet polarizations, current injection from one ferromagnet can, 
due to Andreev reflection, lead to a net drag of current from the second ferromagnet toward the 
superconductor. We present the conditions for the Andreev drag in terms of lead polarizations, 
contact conductances and spin-flip scattering. Remarkably, both equilibrium and nonequilibrium 
zero-frequency current correlations between the ferromagnets become positive even in the presence 
of spin relaxation. 

PACS numbers: 74.45.+C 72.25.-b 72.70.+m 
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I. INTRODUCTION 

At a normal-superconducting interface, Andreev re- 
flection causes a conversion of the quasiparticle charge 
and a flipping of its spin — an electron with spin up, 
incoming towards the superconductor, is reflected as a 
hole with spin down. This leads to exciting conse- 
quences in the transport physics where the focus is put 
on creating, manipulating and detecting spin currents 
and magnetization^^ When the carriers are spin polar- 
ized, the transport properties are strongly sensitive to 
the relative orientation of the ordered moments of fer- 
romagnetic electrodes^ a feature that has received a 
lot of interest due to the relatively long spin dephas- 
ing times observed in metals^ and semiconductors m In 
hybrid systems containing a superconductor, various as- 
pects of the interplay between Andreev reflection and 
ferromagnetism have been addressed theoretically^ and 
exp er iment ally4fl*ii*i^ 

The structure of interest here is shown in Fig.^a): two 
ferromagnetic contacts are connected to a normal con- 
ductor which is in turn connected to a superconductor. 
Of particular interest is the case of crossed Andreev re- 
flection between ferromagnetic leads with antiparallel po- 
larization dominating over crossed normal reflection. In 
such a geometry, with a voltage applied to one of the 
leads while the other lead and the superconductor are 
grounded, the injected current from the biased ferro- 
magnetic lead effectively drags along a current from the 
grounded lead. This effect is easy to understand since the 
transfer of Cooper pairs into the superconductor requires 
a pair of electrons with opposite spin. If the injected 
current consists only of, say, spin-up carriers, the sys- 
tem supplies the spin-down carriers from the oppositely 
polarized ferromagnetic lead even though no voltage is 
applied at this contact. 

The physics of the Andreev drag effect was already 
addressed by Jedema and van Wees et alm& for a multi- 
terminal metallic diffusive conductor. These authors 
pointed to negative four-terminal resistances as a conse- 
quence of the Andreev drag effect. In contrast, much of 
the recent theoretical discussion has focused on a geom- 



etry where two needle-shaped ferromagnetic leads with 
antiparallel polarization are directly coupled to a super- 
conductor at two spatially separated points^ In this 
case, even in the case of 100% spin polarization, the 
crossed Andreev reflection probability is strongly sup- 
pressed in realistic geometries, 1 - making the experimen- 
tal detection of the drag effect difficult. In addition, 
the direct connection of the ferromagnetic electrodes to 
the superconductor might lead to unwanted modifica- 
tions of the local properties of the ferromagneliii and the 
superconductorfli suppressing the drag effect. 

In this work, we overcome the above difficulties by 
inserting a normal multichannel^ mesoscopic conduc- 
tor with generic elastic scattering between the ferromag- 
netic electrodes and the superconductor [see the sketch 
of Fig. ma)]. We include arbitrary interface conduc- 
tances and spin-flip relaxation (see below). The ad- 
vantage of this scheme is twofold: (i) since the super- 
conductor is only contacted at one point, the drag ef- 
fect develops to the fullest extent possible, being of the 
order of the conductances of the contacts between the 
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FIG. 1: (a) Generic ferromagnetic-normal-superconducting 
structure: the process that gives rise to the Andreev drag 
effect is indicated, (b) Sketch of the equivalent circuit model 
(see text). 
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normal conductor and the ferromagnetic and supercon- 
ducting electrodes; and (ii) the sources of the competing 
phenomena — superconductivity and ferromagnetism — 
are spatially separated by the normal conductor. 

We find that the Andreev drag effect is observable 
when the conductance of the normal-superconducting 
(NS) interface dominates over the ferromagnetic-normal 
(FN) contact conductances. The drag effect is reinforced 
with increasing lead polarization and survives for a con- 
siderable amount of spin-flip scattering. Furthermore, we 
demonstrate that the crossed Andreev reflections have a 
profound influence on the zero-frequency current-current 
cross-correlations. We observe that both the thermal 
noise and the shot noise power measured between the 
two ferromagnetic leads yield positive values. 
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II. MODEL 

We analyze a hybrid system (shown in Fig.^l consist- 
ing of a normal conductor connected to two ferromag- 
netic (Fl and F2) and one superconducting (S) reser- 
voirs. For simplicity, the direction of the magnetization 
of the two ferromagnets is assumed to be collinear and 
we consider identical FN couplings. Any type of con- 
tacts, e.g. diffusive-, ballistic- or tunnel contacts, can be 
treated. In what follows, we deal with point contacts at 
the FN (NS) interface with N (M) transversal modes and 
mode independent transparency T (Ts)- We consider the 
case with dimensionless conductances (in units of 2e 2 /h) 
much larger than unity: g = 7VT,<?s = MR 3> 1, where 
R = r|/(2 — Ts) 2 is the Andreev reflection probabil- 
ity. Under these conditions, weak localization as well as 
Coulomb blockade effects can be safely neglected. More- 
over, the dwell time of the normal conductor, T4, is 
assumed to be shorter than typical inelastic scattering 
times. 

Throughout the paper, we consider the case when 
the proximity effect inside the normal conductor is sup- 
pressed (but the Andreev reflection at the NS interface 
is taken into account). A strong proximity effect would 
quench the drag effect as the normal conductor would 
behave effectively as a superconductor mL This can be 
avoided either with a magnetic field (e.g., a stray field 
from the ferromagnets or an externally applied field) or 
with a characteristic quasiparticle energy (bias eV or 
temperature ksT) much larger than the inverse dwell 
time, eV, ksT 3> h/rd- Nevertheless, to avoid single par- 
ticle transport, both eV and ksT should be smaller than 
the gap of the superconductor. 

The resistance of the system is completely dominated 
by the resistances of the FN and NS contacts. More- 
over, we assume that the physical properties within the 
conductor are isotropic due to diffusive or chaotic scat- 
tering. We can thus consider the normal conductor to be 
effectively zero dimensional. This restricts our model to 
systems with a spin-flip length much larger than the lat- 
eral dimensions of the system, i.e., there are no spatial 



FIG. 2: Off-diagonal conductance G21 = {h/2e 2 )T 2 /Vi of the 
current in ferromagnetic lead F2 when a voltage bias is applied 
to lead Fl as a function of the contact conductances. The 
magnetic moments of the reservoirs are aligned parallel in (a) 
with spin-flip rate 7 = and (b) with polarization 77 = 1. (c) 
and (d) same as (a) and (b) in the antiparallel orientation. 



variations of accumulated spins inside the dot, though 
no assumptions are made about the relation between 
the spin- flip time r s f and r<j. We stress that all the 
above constraints form a situation that is experimentally 
accessible with present techniques. Vertical structures 
like metallic thin-film heterostructures, epitaxial layers of 
half- metallic chromium dioxide or diluted magnetic semi- 
conductors are equally suitable for the observation of the 
effects discussed here. In fact, a uniform spin accumu- 
lation was very recently achieved in a normal metallic 
island coupled to magnetic leads i& 

General theoretical formulations of semiclassical spin 
transport, based on quasiclassical Green's functions, have 
been developed by Brataas, Nazarov and Bauer for 
normal-ferromagnetic systems 1 - and by Belzig et al. for 
normal-ferromagnetic-superconducting systems^ Here, 
we follow the approach of Samuelsson and Biittiker— ex- 
tended to spin-dependent transport. This simple theory 
is equivalent to Refs. I9IIT9I for the geometry in Fig^a). 
It is directly formulated in terms of distribution func- 
tions and allows us to treat the current and the current 
correlations within a unified framework. Because the 
proximity effect is suppressed, the electrons and holes 
are independent^ and the system is mapped onto the 
equivalent circuit of Fig. ^b) in which the electron and 
hole subsystems (with opposite spins) are connected via 
Andreev reflection. Spin relaxation in the normal con- 
ductor gives rise to a spin-flip current between the spin 
subsystems with the same quasiparticle charge. As a re- 
sult, we can describe the conductor with four nodes [two 
nodes (t and j) for each quasiparticle type (electrons 
and holes)] connected to each other and to the ferromag- 
netic reservoirs via spin dependent resistances. We model 
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the spin dependent conductances of the FN contacts as 
9a = -9(1 ^ 7 ?ct)/2, where rj a is the spin polarization of 
the lead a (a — {1, 2}): rj a = {v\ - v a )/{v a + v^), where 
i/£ is the density of states of the ferromagnetic lead close 
to the Fermi energy (a = {1,1}). 

Conservation of time-averaged (spin-dependent) cur- 
rent densities must be fulfilled at each node for every 
energy E. Currents flowing into the normal conductor 
are defined positive such that the spectral currents obey 

J2 ^1 Ti S + Ht,e = > T - HT,e = . C 1 ) 

a a 

where q = (e, h) labels the carrier (electron or hole) and 
— (+) is for q = e (q = h). The currents through the FN 
and NS contacts are: 

? aq = g a a {U - r q ), r s = 9s(J!-ft)/2, (2) 

where f aq (E) = {1 + cxp ([E T eV^/feT)} -1 is the equi- 
librium distribution function at reservoir a. Note that 
f aq is spin independent as spin relaxation mechanisms 
are usually highly efficient in ferromagnets. In Eq. J5J, 
fq(E) is the nonequilibrium distribution function of the 
normal conductor to be calculated from Eq. QJ. The 
spin-flip current i^ q — j(f^ — /J) equilibrates the spin 
distributions. Here, 7 = hfo/T s t (vq is the density of 
states of the conductor at the Fermi energy and r s f the 
spin- flip time) 

III. CONDUCTANCE MATRIX 

The total current flowing out of reservoir a is 

where T£ = (e//i) J dEi^ (E), from which we determine 
the conductance matrix. The conductance G a j3 relates 
the current X a with a voltage change at lead f3. In the 
case V\ > and V% = Vs = 0, the interesting conduc- 
tance is G 2 i = {h/2e 2 )(I 2 /V 1 ): 

q = g 2 [g(vi 2 + m 2 - 2) - 4 7 - 2 mmgs \ 

21 .9 2 [4 - (vi + V2) 2 } + 8755 + 4 5 (2 7 + gs) ' 1 ' 

In a normal system, G21 is always negatived On the 
contrary, in our system we can have G21 > 0, which is 
the manifestation of the Andreev drag effecti^ Figure [21 
shows the behavior of G21 as a function of 2g/gs- In 
the parallel (P) case, r/i = 772 = r\ and G21 is nega- 
tive [Fig. |2f a)]. The maximum (absolute) value of G21 
is attained at 77 = 1 and decreases with increasing values 
of 7 [Fig. |21[b)] . In the antiparallel (AP) configuration 
(r/i = — 7/2 = 77), a reversal of the current at F2 occurs 
for low values of g/gs and 7 [see Pig- EJc) and Fig.^d)]. 
Note that the drag effect can result already for a small 
value of 77. In fact, crossed Andreev reflections start to 
dominate when if > (g/gs) /(l + 9/ 9s) a t 7 = 0. In 



particular, for the fully polarized case one always finds 
G21 > regardless of the ratio g/gs- When gs ^> g, we 
find that G21 > holds for any value of 77 > 0, in agree- 
ment with Ref. Increasing 7, the spin relaxation pro- 
cesses randomize the spin distributions and eventually 
lead to G21 < 0, [the curves tend to the 77 = case of 
Fig. Gib)]. The condition for the drag effect is: 

27<^ 2 3s-.9(l-^ 2 ). (5) 

This is a central result of our work. 

Interestingly, the Andreev drag effect implies that even 
the equilibrium Nyquist- Johnson noise yields positive 
current cross-correlations P21 = 4fcsTG2i in contrast to 
normal conductors. Nevertheless, energy dissipation is 
ensured as the eigenvalues of the conductance matrix are 
semipositive. 

Another striking result is that for 77 = 1 there are two 
linearly independent voltage configurations that yield 
zero eigenvalue of the conductance matrix. One is the 
gauge solution (V\ = V2 = Vs). The nontrivial arrange- 
ment is V\ = — V2 = V and Vs = 0. Here, Fl injects elec- 
trons whereas an equivalent flow of incident holes stems 
from F2. Complete Andreev reflection yields then zero 
net current. 



IV. MAGNETIZATION 

Out of equilibrium for V\ > and V2 = Vs — 0, an 
accumulation of majority spins forms within the conduc- 
tor. Before applying the voltage, the conductor has zero 
magnetization. In the AP case, the polarization excess is 
given by 

M _ JdEE q (P q -n) _ gr, 

JdEJ2 q W,+U) 2 7 + .g + ffs ' 

It approaches 77 for g 3> gs , 7 and vanishes for 7 ^> g or 
9S > 9- 

V. CROSS CORRELATIONS 

To obtain additional insight into the drag effect, we in- 
vestigate the shot noise of the current. Of particular in- 
terest are the correlations between currents flowing in the 
two ferromagnetic reservoirs. In purely normal systems, 
the current cross-correlations are manifestly negative. 23 
However, in hybrid superconducting structures it was 
shown that Andreev reflections can cause a reversal of 
the sign of the cross-correlations»24 The positive correla- 
tions in few-mode conductors were also shown to be en- 
hanced by ferromagnetic leads^ Here, we focus on the 
cross correlations in the AP configuration, which is the 
most interesting case. 

The contacts act as generators of fluctuations 6i° q and 
Sig. As a result, the spin-dependent distribution func- 
tion of the normal conductor fluctuates itself. The total 
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FIG. 3: Cross correlation P12 for T = 1 and R = 1 versus 
the ratio of the conductances of the point contacts in the 
antiparallel case. In (a) 7 = whereas in (b) 77 = 1. Insets: 
"Phase diagram" of the occurrence of the Andreev drag effect 
(below the dashed line) and positive P12 (below the solid line) 
for (a) 7 = and (b) 77 = 1. 

fluctuating current is 
^ aq = K q -9 a Jf^ ^ s = Ws+9s{5f s e -8f° h ). (7) 

We also take into account the fluctuations Sitf q due to 
spin relaxation, giving the fluctuating current 

A lLhq = Si Lhq + 1 (5ft-5p q ). (8) 

The next step is to determine the fluctuations of the 
distribution functions of the conductor, 5f q . In the zero 
frequency limit, the fluctuating currents are conserved 
(=F for q = e,h): 

^A4,tA4+Ah u = 0, (9a) 

a 

^A4tA4-Ah u = 0. (9b) 

a 

The fluctuations of the charge current at contact a are 
found from Ai a = Ai } ae + Ai^ e - Ai ] ah - Ai l ah . 
The noise power of the current cross-correlations is 

P21 = 2 J dt (AX 2 (t)Ari(0)) . (10) 

The contacts emit fluctuations independently so that 

K,(B,t)^;,,(£',o)^<, (11) 

with 5 = 8 aa ,8 qq ,5 aa ,5{E - E')S(t - t'). Here, S° q is the 
fluctuation poweri22iS 

S° aq = eg°[fZ q + r q - 2f°r aq - T{f° q - f-f] ,(12a) 
S° s = eg s [t + fl - 2 fit - R(fe ~ flYM^ (12b) 
S iU = e 7 [/j + fl - 2/J/Jr] , (12c) 



where we have assumed that the statistics of spin-flip 
events is Poissoniani2& 

The above system of equations is complete and can be 
solved for any set of applied voltages, polarization align- 
ments and spin-flip rates. We calculate P21 for the bias 
configuration given by V\ > 0, Vi = Vs = in the AP 
case at ksT = 0. We present only the analytical result 
in limiting cases and show the general result in Fig. [3] 
In Fig. 12 a), we plot P21 as a function of the conduc- 
tance ratio 2g/g$ for various polarizations 77, 7 = and 
r = R = 1. (Qualitatively, the results are independent 
of r and R.) For 77 = 0, the cross-correlations are nega- 
tive. As 77 is turned on, P21 can become positive. With 
increasing g/gs, P21 exhibits a crossover from negative 
to positive values at 2g/gs ~ 1. Indeed, for dominating 
NS conductance (small g/gs) we find 

P21 - PoV 2 [r + 2 - ?7 2 (3r - 2)]/T6 + 0(g/g s ) , (13) 

which is manifestly positive for 77 > (po = &gVie / h) . 
In this case, the electron and hole distributions within 
the conductor are almost identical. 

On the other hand, for dominating FN conductance 
(small gs/g) the cross-correlations become 

P21 =p a (l-rf)[T-2 + 7f(3T-2)}/8 + 0(g s /g), (14) 

which are negative for rj < 1. Notice that for 77 = 
this expression yields the shot noise for a chaotic dot 
coupled to two normal leads. In the fully polarized case, 
the cross-correlations are always positive independently 
of g/gs- 

For illustrative purposes we show in the inset of 
Fig. the "phase diagram" (rj,g/gs). The solid line 
marks the crossover between positive and negative cross- 
correlations. With increasing 77, the crossover point shifts 
toward larger values of 2g/gs- For comparison, we also 
plot a dashed curve that marks the transition to the An- 
dreev drag effect. Notice that the conditions P21 > and 
G21 > are closely related, but not identical. This is the 
case since the scattering processes in general contribute 
in different ways for the current and noised 

A finite source of spin relaxation leads to the suppres- 
sion of positive P21. To see this, we show in Fig.[3{b) the 
cross-correlations for the fully polarized case for different 
values of 7. For nonzero values of 7, P21 becomes neg- 
ative with increasing 2g/gs- Eventually, when 7 — > 00, 
we recover the unpolarized case (77 = 0) of Fig. 02a). 

VI. CONCLUSIONS 

We have investigated the transport properties of a nor- 
mal conductor attached to ferromagnetic leads and cou- 
pled to a superconductor in the subgap regime. The pro- 
posed geometry, which includes generic elastic and spin- 
flip scattering, is shown to exhibit a pronounced Andreev 
drag effect. We have demonstrated that the equilibrium 
and transport current correlations are strongly sensitive 
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to crossed Andreev reflections. The theory presented 
here applies to a wide range of geometries and should 
therefore advance experimental observation of the An- 
dreev drag effect. 

Note added in proof. Recently, we became aware of the 
work of Lambert et alm^ who compute numerically the 
conductances for a metallic diffusive geometry similar to 
Fig. 1. They also conclude that a normal conductor in- 
serted between the ferromagnetic and the superconduct- 
ing contacts leads to an enhancement of the Andreev drag 
effect. In contrast to our work, the role of spin- flip, the 
effect of differing lead polarizations, and the shot noise 



are not investigated. 
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